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Abstract 

MACHe3 (MAtrix of Cells of superfluid '^He) is a project of a new detector for direct 
Dark Matter search, using superfluid ^He as a sensitive medium. This paper presents 
a phenomenological study done with the DarkSUSY code, in order to investigate 
the discovery potential of this project of detector, as well as its complementarity 
with existing and planned devices. 

Key words: Dark Matter, Super symmetry, Superfluid Helium-3, Bolometer. 
PACS : 95.35; 67.57; 07.57.K; 11.30.P 



1 Introduction 



A substantial body of astrophysical evidence supports the existence of non- 
baryonic Dark Matter (DM) in the halo of our galaxy, in particular in the form 
of new, yet undiscovered, weakly interactive massive particles (WIMPs) [1]. 
One of the leading candidates is the neutralino predicted by the supersym- 
metric extensions of the Standard Model of particle physics. 
Following early experimental works [2], a superfluid ^He detector has recently 
been proposed [3] for direct Dark Matter search. Monte Carlo simulations have 
shown that a high granularity detector, a matrix of superfluid ^He cells, would 
allow to reach a high rejection factor against background events, leading to a 
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low false event rate. The purpose of this paper is to present a full estimation of 
the neutralino (x) event rate, within the framework of the Minimal Supersym- 
metric Standard Model (MSSM), in order to compare with background event 
rate obtained by Monte Carlo simulation [3,4]. Finally, the complementarity 
with existing devices, both for direct and indirect detection, will be shown. 

1.1 Experimental device 

The elementary component of MACHeS is the superfluid ^He cell [5,6]. It is 
a small copper cubic box (V ~ 125 mm^) filled with superfluid ^He in the 
B-phase. This ultra low temperature device (T ~ 100 /iK) presents a low de- 
tection threshold {Eth — 1 keV). An experimental test of such a prototype cell 
has been done at CRTBT in June 2001. Preliminary results [7,8] show that 
a threshold value down to ~ 1 keV has been achieved, and that a stability of 
the order of one week, at T ~ 100 /iK, has been obtained. 
The final version of the detector will be a matrix of 1000 cells of 125 cm^ 
each. The idea is to take advantage both on the energy loss measurement and 
the correlation among the cells to discriminate neutralino events from those 
of background (neutrons, 7-rays and muons). The design of the matrix has 
been optimized, with a Monte Carlo simulation [3]. In the preferred configu- 
ration, a 10 kg detector, the false event rate has been shown to be as small 
as ~ 10^^ day~^ for neutron event^ and ~ 10~^day~^ for muon events [3]. 
Background from 7-ray events needs to be taken into account. Energy loss 
measurement and correlation among the cells within a 10 kg detector allows 
to obtain a rejection up to 99.8% for ~ 2 MeV 7-rays [3]. Additionnal internal 
tag on 7-rays may be obtained using a new matrix configuration in which two 
neigbouring cells share a common copper wall, thus greatly improving corre- 
lation factor while reducing the amount of copper used. Monte Carlo studies 
are under way [8]. Furthermore, 7-ray contamination is to be estimated in the 
forthcoming months for a multicellular prototype matrix. 
In order to compare these false event rates with the expected neutralino rate, 
only muons and neutrons will be taken into account in the following. We shall 
recall that neutrons are usually considered as the ultimate background noise 
for this type of search, as they interact a priori like WIMPs. 

1.2 ^He as a sensitive medium for direct DM search 

Several properties of ^He make this nucleus a promising candidate for a sen- 
sitive medium for direct DM search. 

^ Fast neutron contribution, from interaction of muons in the rock, is expected to 
be negligible as the enregy release in the cell will be much greater than 6kev. 



2 



a) Concerning background rejection, and as outlined in [3,9], the neutron 
capture process offers the possibihty to discriminate neutron and x event, 
when considering a 10 leg granular detector. Compton cross-section being small 
(cr < Ibarn), the interactions with 7-rays will be minimized. Eventually, as 
explained in [6], superffuid ^He is produced with an extremely high purity, the 
only solute being ^He in a negligible fraction. Consequently, no contamination 
from radioactive materials is expected in the sensitive medium. Of course, 
natural radioactivity from external materials (Cu, ...) has to be taken into 
account, by a careful selection of these materials. 

b) Concerning neutralino detection the advantage is twofold. First, the max- 
imum recoil energy does only slightly depend on the neutralino mass, due to 
the fact that the target nucleus (m = 2.81GeV/c^) is much lighter than the 
incoming x i^x — 32GeV/c^), considering latest results from collider experi- 
ments [10]. As a matter of fact, the recoil energy range needs to be studied only 
below 6keV, see [3,7]. Second, ^He being a 1/2 spin nucleus, an ^He detector 
will be sensitive mainly to axial interaction, making this device complemen- 
tary to existing ones, as shown below. In fact, the axial interaction is largely 
dominant (up to three orders of magnitude) in all the SUSY region associated 
with a substantial elastic cross-section [7]. 

1.3 Theoretical framework 

This phenomenological study has been done with the DarkSUSY codeQ [11], 
within the framework of the phenomenological Supersymmetric Model, namely 
with the following free parameters : 

/i, M2, tan/3, ruA, rriQ and Af, t (1) 

with /J, (M2) the Higgsino (Gaugino) mass parameter, tan (3 the ratio of Higgs 
vacuum expectation values, tua the CP-odd Higgs boson mass, mo the com- 
mon scalar mass and Ah t the soft trilinear coupling parameters^. 
Apart from A^ t chosen at fixed zero value, as their influence is expected to 
be negligible, all the parameters have been scanned on a large range, with a 
variable number of steps (tab. 1). This scan of the free parameters corresponds 
to a total number of supersymmetric (SUSY) models of the order of 2 x 10^. 
We shall suppose all through this work the neutralino (x), the lightest su- 
persymmetric particle, as the particle making up the bulk of galactic cold 
DM. Each SUSY model is then checked not to be excluded by collider exper- 
iments [10], including 6 ^ 57 limit. 

^ The version used is 3.14.01, with correction of some minor bugs. 
For a good introduction to MSSM models, we refer the reader to [12]. 
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The nest step, the evaluation of the rehc density, is the key point of any Dark 
Matter calculations. Given the number of free parameters (5 in our case), the 
allowed SUSY parameter space may be extremely large, leading to a x relic 
density ranging on up to five orders of magnitude [7]. In order to exclude 
SUSY models giving a x relic density too far away from the estimated matter 
density in the Universe [13] [Qu — 0.3), only models with in the following 
range are considered : 

0.025 < fi^ho^ < 1 (2) 

where ho = (0.71 ± 0.07) Xqqj is the normalized Hubble expansion rate [10]. 
The lower limit comes from the condition that the neutralino relic density has 
to be at least greater that the baryonic density, and the upper limit is a con- 
servative limit so that x do not give a density greater than the UniverseQ. The 
SUSY model is thus checked to really provide a good non-baryonic Dark Mat- 
ter candidate. We follow [14] in the choice of the " cosmologically interesting" 
range of fl^. It should be noticed that this loose selection allows to exclude a 
large number of models in our SUSY scan. For further details concerning the 
calculation of the x relic density, we refer the reader to [11]. 

As the detection on earth is concerned, a galactic halo model has to be con- 
sidered. In the following, standard parameters have been used, in a spherical 
isothermal halo distribution, with a local density (po) and an average veloc- 
ity (fo), with the following values : 

Po = 0.3 GeV/c^ cm~^ and vq = 220 km s"^ (3) 

These parameters are widely used for dark matter detection computations, see 
[15] for instance. No clumpy galactic dark matter structures [16] are considered 
hereafter, since the effect is expected to be small both for direct detection and 
neutrino telescopes, in which the signal depends on the local halo density, as 
emphasized in [17]. 

2 Spin dependent cross-section and event rate 

As early recognized by Goodman and Witten [18], the interaction between 
neutralino and quarks may be either spin dependent or spin independent, in- 
volving different Feynman diagrams [1]. In the general WIMP case, the allowed 
interactions are : vectorial, axial and scalar. The neutralino being a Majorana 
fermion, the vectorial interaction vanishes, leaving two classes of interaction : 
scalar and axial, the first one being spin independent, the second spin depen- 
dent and obviously requiring a non-zero spin nucleus. ^He being a light 1/2 

^ It can be noticed that selecting on ri^ho^ allows for a slightly looser selection. 
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spin nucleus, a medium made of such nuclei will be sensitive only to axial 

interaction, as shown in [7]. 

Using the DarkSUSY code, the x-^He spin dependent cross-section has been 
evaluated. The calculation of the %- quark elastic scattering amphtude is done 
at the tree- level, via an exchange of squark or Z^. The amplitude on nu- 
cleon (ap/n) is then cvahiatcd by adding the contribution of each quark (a^j), 
weighted by the quark contents of the nucleon : 

A(p/n) . 



where is the quark contents of the nucleon [19]. 

The axial cross-section on ^He is then given by : 

32 (J + 1) 

aspinCHe) = — Gl,ml — (op <Sp> + an<Sn >f (5) 



where < Spjn > is the spin contents of the ^He nucleus (< Sp >— —0.05 and 
< Sn >— 0.49), rur is the reduced mass and J the ground state angular mo- 
mentum of the ^He nucleus. 

Figure 1 presents the cross-section on ^He as a function of the x mass. It 
can be seen that, for SUSY models not excluded by coUider experiments and 
giving a relic density within the range of interest, a cross-section as high as 
~10~^pb can be obtained for ~60GeV/c^ neutralino. 

Using the cross-section, the x event rate (R) has been evaluated for a 10 kg 
^He matrix, for comparison with the background rate previously evaluated [3] . 

with pQ the local halo density and Vq the average neutralino velocity (eq. 3). 
A large number of models are giving a rate higher than the estimated false 

event rate induced by neutrons ( ~ 0.1 day~^), or above the estimated muon 
background (~ 10^^ day^^). It could thus be concluded that a high granularity 
^He detector would present a sensitivity to a large part of the SUSY region. 

In the following, the ji background level (lO""^ day~^) is considered to be the 
lowest reachable limit for MACHe3 and it is thus taken as the reference value. 
Any model giving a rate greater than this value is considered hereafter to be 
visible. This will be used for the comparison with other DM search strategies. 
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3 Complementarity with existing devices 

We present a study of the discovery potential of various detection strate- 
gies that may be correlated with MACHeS. It is indeed worth understanding 
whether planned projects would be sensitive to different SUSY regions. In each 
case, we will take into account background estimations, or projected limits, 
in order to investigate the different discovery potentials. Choosing conserva- 
tive values will allow to study the phenomenology of various searches without 
bias. It should be noticed that the results are model-dependent, both SUSY 
and galactic halo ones. However, standard galactic parameters have been con- 
sidered and the large SUSY scan allows for an exhaustive study within phe- 
nomenological MSSM model. 

3.1 Complementarity with scalar detection 

Firstly, we studied the possible correlation between the two types of direct 
detection : axial and scalar. In the latter case, many detectors are already 
running, or planned to start measurements in the near future. At first sight, 
these two direct searches should be largely independent as they involve dif- 
ferent processes, in particular the Feynman diagrams at the tree level are 
different. 

Within the framework of the study described above, the proton scalar cross- 
section has been evaluated at the tree levelQ. Figure 2 presents the scalar 
cross-section on proton as a function of the x mass, for all models of the scan 
not excluded by collider experiments and associated with a relic density within 
the range of interest (eq. 2). It can be seen that many SUSY models lie below 
the projected limits of future scalar detectors (CDMS [21], CRESST [23]), 
while giving an event rate above the value taken as the lowest reachable limit 
for MACHeS. 

As an illustration of this complementarity, we present on fig. 3 the regions 
potentially covered by direct scalar search and by MACHeS, in a given part 
of the SUSY parameter space (see fig. 3 for details). The reference limits have 
be chosen as : the CDMS projected limit [15] in one case and the 10~^ day~^ 
limit for MACHeS, previously discussed. 

It can be noticed that only axial direct search may be sensitive to the /i < 
region, in this given part of the SUSY parameter space, due to the fact that 
scalar interaction vanishes. The cancellation arises in the t channel (Higgs 
exchange) which is dominant, with u and d quark contributions of opposite 
sign and absolute values of the same order [24]. Similar cancellation does arise 

^ Higher order diagrams, involving gluon loops [1], are not taken into account in 
the current version of DarkSUSY. 
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for different value of tan /9, in the case of axial interaction. It can be concluded 
that these two detection strategies are sensitive, as expected, to different SUSY 
regions, thus highlighting their complementarity. 



3.2 Complementarity with indirect detection 

Amongst various kind of indirect searches (p, D, e"*", 7, u), the detection ex- 
pected in neutrino telescopes is the only one likely to be correlated with di- 
rect detection, as it involves both elastic scattering and annihilation. The neu- 
tralino capture may occur either in the Sun or in the Earth, through successive 
elastic scatterings. As far as comparison with direct detection is concerned, 
the key point is that the composition of the Sun and the Earth are different. 
While the Earth is mainly made of even-even nuclei (O, Si, Mg, Fe), the Sun is 
mainly made of H and He. So one expects the signal coming from annihilation 
in the centre of the Earth to be only correlated with the signal expected from 
scalar direct detection. On the contrary, the Sun being made of both even 
and odd nuclei, a very light correlation is expected with direct detection, both 
scalar and axial. 

Within the framework of the study described above, we tried to find out 
whether the MACHeS project would present a complementarity with the 
project of u telescopes, such as IceCube [26] or Antares [27]. These projects 
aimed at building a very large detector, of size of the order of Ikm^. The 
instrumentation will thus be loose, giving a rather high detection threshold. 
We follow [25] in choosing a reference threshold of 25 GeV. 

Figure 4 shows the result of the calculation for u telescopes and MACHeS. The 
expected n flux (km^^year^^) is shown against the rate in MACHeS (day^^). 
The background usually taken as ultimate for 1/ telescopes, coming from the 
interaction of cosmic rays in the sun's corona [28], is indicated at a level of 
~ 10 km~^ year~^. We consider this value as the most favourable limit for these 
projects of km'^ detectors. In the same way, the background noise from /i in 
MACHeS is taken as the lowest reachable level, as discussed above (sec. 2). 
Any model above these limits are said hereafter to be visible for this type of 
detection. It can be concluded from fig. 4 that there is a complementarity be- 
tween indirect u detection and axial direct detection, and in particular in the 
case of the MACHeS detector. Many models may be visible in only one kind 
of detection strategy, high mass neutralinos being seen, for instance, only in 
u telescopes. 

As for the comparison between scalar and axial detection, we present on fig. 5 
the SUSY regions potentially observable for these two types of detection strate- 
gies. As expected, neutrino telescopes are sensitive to higher values of M2. For 
this given region of the SUSY parameter space, it can be concluded that the 
combination of these two signals may close the gap between collider excluded 
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region and cosmological bounded regions, up to high M2 values. 
4 Conclusion 

It has been shown that a 10 kg high granularity ^He detector (MACHeS) would 
allow to obtain, in many SUSY models, a x event rate higher than the esti- 
mated (neutrons and muons) background. MACHeS would thus potentially 
allow to reach a large part of the SUSY region, not excluded by current col- 
lider limits and for which the neutralino relic density lies within the range 
of interest. Furthermore, it has been shown that this project of new detector 
would be sensitive to SUSY regions not covered by future or ongoing DM 
search detectors, both for direct and indirect detection, thus highlighting the 
complementarity of MACHeS with existing or planned devices. 
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Table 1 

Scan of the SUSY parameters used for the study. The total number of models is of 
the order of ~ 2 x 10^. 
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-1 0.025<Qh^<1 DarkSUSY 




Fig. 1. Axial cross-section on ^He (pb) as a function of the neutralino mass (Gev/c'^), 
for all models of the scan not excluded neither by colliders nor by cosmological 
constraints (eq. 2). The rates corresponding to the Monte Carlo evaluation [3] of 
neutron and muon background are indicated for reference. 
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Neutralino mass (GeV/c^) 



Fig. 2. Scalar cross-section (on proton) as a function of the x mass, for all models of 
the scan not excluded neither by colliders nor by cosmological constraints. Exclusion 
limits from the Edelweiss [20] and CDMS [21] experiments are shown, as well as the 
3a "DAMA region" [22]. Dotted lines indicate projected limits [15] from CRESST 
[23] and CDMS. Dark points indicate SUSY models giving a x rate in MACHeS 
lower than the estimated background level, while light points are giving a rate 
higher than 10~^ day~^. The two regions are overlapping as the signal in MACHeS 
is not correlated with scalar cross-section. 
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Fig. 3. Sensitivity regions for scalar search and MACHeS, in the (M2,/x) plane. The 
three figures have been obtained with tan (] = 10, mg — 1.6 TeV/c^ and three values 
of Ma- The central white region is excluded by collider limits, while the two other 
white regions are associated with a relic density outside the preferred range. The 
color code indicates the region potentially covered by MACHeS (with a 10~^ day~^ 
background value), by scalar detection (for which CDMS projected limit is taken 
as ultimate limit), by both strategies, or none of them. 
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Fig. 4. Muon flux in neutrino telescopes (with a 25 GeV threshold) against event 
rate in MACHe3. For references, the estimated backgrounds are shown, in one case 
from energetic neutrinos created in the sun's corona [28], and in the other case the 
background from muon and neutrino events in MACHeS. 
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Fig. 5. See fig 3. In this case the color code indicates the region potentially covered 
by MACHeS (with a 10~^ day~^ background value), by neutrino telescopes detection 
(for which a 10 km~^ year~^ background value [28] is considered), by both strategies, 
or none of them. 
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